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■ Abstract. The cosmological relic density of the lightest supersymmet- 
I I ric particle of the minimal super symmetric standard model is calculated 

under the assumption of gauge and Yukawa coupling unification. We em- 
ploy radiative electroweak breaking with universal boundary conditions 
from gravity-mediated supersymmetry breaking. Further constraints are 
J> I imposed by the experimental bounds on the 6-quark mass and the BR{b — > 

0^ ■ 57). We find that coannihilation of the lightest supersymmetric particle, 

which turns out to be an almost pure bino, with the next-to-lightest su- 

■ persymmetric particle (the lightest stau) is crucial for reducing its relic 
O . density to an acceptable level. 
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1. Introduction 

D , In recent years the consideration of exotic dark matter has become necessary in 

' order to close the Universe In the currently favored supersymmetric (SUSY) 

extensions of the standard model, the most natural WIMP and candidate for 
CDM is the LSP, i.e. the lightest supersymmetric particle. In the most favored 
scenarios the LSP is the lightest neutralino, which can be simply described as 
d I a Majorana fermion, a linear combination of the neutral components of the 

gauginos and Higgsinos |lO|. Its stability is guaranteed by imposing R-parity 
conservation, which implies that the LSP's can disappear only by annihilating 
in pairs. 

The simplest and most restrictive version of MSSM with gauge coupling uni- 
fication is based on the assumption of radiative electroweak symmetry breaking 
with universal boundary conditions from gravity-mediated soft SUSY breaking. 
An interesting question is whether this scheme is compatible with exact "asymp- 
totic" unification of the three third family Yukawa couplings. A positive answer 
to this question would be very desirable since it would lead to a simple and 
highly predictive theory. 

In this presentation we summarize the results of Refs. [^, ^, aimed to answer 
the question of whether the simple version of the MSSM with Yukawa unification 
at the GUT scale, can be compatible with the most restrictive phenomenological 
constrains, which are the correct predictions for 6-quark mass and BR(b 57), 
and satisfy the requirement that the relic abundance ^lsp of the lightest 
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super symmetric particle (LSP) in the universe does not exceed the upper limit on 
the cold dark matter (CDM) abundance implied by cosmological considerations. 

2. Input parameters on the MSSM with Yukawa Unification 

We consider the MSSM embedded in some general supersymmetric GUT based 
on a gauge group such as 50(10) or Eq (where all the particles of one family 
belong to a single representation) with the additional requirement that the top, 
bottom and tau Yukawa couplings unify at the GUT scale Mqut- Ignoring 
the Yukawa couplings of the first and second generation, the effective superpo- 
tential below Mqut is 

W = eiji-htHiQif + h^HlQib^ + KHILIt^ + ^lH\Hi) , (1) 

where Qs = (t, 6) and L3 = {ur,T) are the quark and lepton SU{2)l doublet left 
handed superfields of the third generation and f^, and the corresponding 
SU{2)l singlets. Also, Hi, H2 are the electroweak higgs superfields and £12 = 
+1. The gravity- mediated soft supersymmetry breaking terms in the scalar 
potential are given by 

a,b 

(eiji-AthtHiQii" + AbhHiQib'' + ArhrH{L{f^ + BfiHlHl) + h.c.) , (2) 

where the (pa 's are the (complex) scalar fields and tildes denote superpartners. 
The gaugino mass terms in the Lagrangian are 

-,3 8 

- (MiBB + M2 J2 WrWr + M3 ^ gaga + h.c.) , (3) 
r=l a=l 

where B, Wr and ga are the bino, winos and gluinos respectively. 'Asymptotic' 
Yukawa coupling unification implies 

htiMcuT) = hbiMcuT) = hriMcur) = h . (4) 

Based on iV = 1 supergravity, we take universal soft supersymmetry breaking 
terms at Mgut, i-e., a common mass for the scalar fields mo, a common trilinear 
scalar coupling Aq and Bq = Aq — niQ. Also, a common gaugino mass M1/2 is 
assumed at Mqut- 

Our effective theory below Mqut depends on the parameters (/^o = fJ-iMGUT)) 
mo, Ml/2, ^0, /"o, OG, Mgut, ho, tan/3 . 

ac and Mgut 

are evaluated consistently with the experimental values of o-efyi, o-s 
and sin^ ^vK at "mz- We integrate numerically the renormalization group 
equations (RGEs) for the MSSM at two loops in the gauge and Yukawa 
couplings from Mgut down to a common supersymmetry threshold Ms = 
yjvni_^ rrii^ . From this energy to mz, the RGEs of the nonsupersymmetric 
standard model are used. 
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tan/3 is estimated at the scale Ms using the experimental input mr{m-j 

1.777 GeV. We incorporate the SUSY threshold correction to m.j-{Ms) 



from the approximate formula of Ref.|14]. It is about 8%, for > 0, 
leading to a value of tan /? = 55.4 — 54.5 for rriA = 100 — 700 GeV, while, 
for /i < 0, we find a correction of about —7% and tan/3 = 47.8 — 46.9 in 
the same range of tjia- 

ho is found by fixing the top quark mass at the center of its experimental range, 
^tiiTit) = 166 GeV. The value obtained for mb{mz) after including su- 
persymmetric corrections is somewhat higher than the experimental limit. 

Aq , for simplicity we take = 0. Our results move very little for negative 
values of Aq bigger than about —.5Mi/2, however lower negative values 
and positive values of this parameter tend to increase the SUSY spectrum 
increasing i^iSP h'^- Therefore the limits on 0,lsp imposes lower and 
upper bounds on Aq. 

niQ, Mi/2, fJ'O, As we will discuss later the electroweak symmetry breaking, when 
Yukawa unification is assumed, imposes a relation on the values of mo and 
Ml/2- On the other hand the role of coannihilation make convenient 
to express our results in terms of relative mass splitting Af^ = {rrif^ — 
m^/m^ between the NLSP and LSP. Therefore we trade the GUT values 
of mo. Ml 12 and /io by the pseudoscalar Higgs mass niA and Afj. 

Let us describe with more detail the last item above. Assuming radiative elec- 
troweak symmetry breaking, we can express the values of the parameters ^ (up 
to its sign) and B at Ms in terms of the other input parameters by means of 
the appropriate conditions 



tan2 (3-1 ^rnz^^mip- ^3^^ ^ ^^^^ ^ 



where m//^, m/^j the soft supersymmetry breaking scalar higgs masses. 
When unified Yukawa couplings and a common value for m/^j and are 
assumed GUT , we find both m^^^ and m|^^ negative at Ms- However for cer- 
tain values of mo and is possible to find values the pseudoscalar Higgs, 

vriA = mjj^ + mjj^ + 

beyond a lower bound, considered to be mz in the present work. Furthermore, 
the authors of Ref. |^ found that, for every value of ruA and a fixed value of 
mt{mt), there is a pair of minimal values of mo and M1/2 where the masses of 
the LSP and f2 are equal. This is understood from the dependence of niA on 
mo and Mx/2 given in Ref. |^: 

m\ = aMf/2 ~ Pf^o ~ const. , (6) 

where all the coefficients are positive and a and /3, which depend only on mt{mt), 
are ~ 0.1 (the constant turns out to be numerically close to m^). Equating the 
masses of the LSP and f2 is equivalent to relating mo and Mi/2- Then, for every 
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Figure 1. The values of m^, rriQ, M1/2 and Ms as functions of rriA 
for /i > 0, j4o = and = rny-. These values are affected very little 
by changing the sign of fi. 

vriA, a pair of values of rriQ and Mi 12 is determined. We had included the full 
one-loop radiative corrections to the effective potential as given in the appendix 
E of Ref. 0. 

The values of the LSP, Ms and the corresponding values of mo and M1/2 
are given in Fig.l. 

3. Phenomenological constraints from and 6 — > 57 

A significant problem, which may be faced in trying to reconcile Yukawa unifi- 
cation and universal boundary conditions, is due to the generation of sizeable 
SUSY corrections to the 5-quark mass H, |^]. The sign of these corrections is 
opposite to the sign of the MSSM parameter // (with the conventions of Ref.Q). 
As a consequence, for < 0, the tree-level value of m;,, which is predicted 
from Yukawa unification already near its experimental upper bound, receives 
large positive corrections which drive it well outside the allowed range. How- 
ever, it should be noted that this problem arises in the simplest realization of 
this scheme. In complete models correctly incorporating fermion masses and 
mixing, mf, can receive extra corrections which may make it compatible with 
experiment. So, we do not consider this 6-quark mass problem absolutely fatal 
for the // < case. However in the alternative scenario, with /U > 0, the &-quark 
mass receives negative SUSY corrections and can easily be compatible with data 
in this case. An example of the typical values we find is given in table I. 

This scheme with /x > 0, is severely restricted by the recent experimental 
results on the inclusive decay 6 — > 57 |^. It is well-known that the SUSY 
corrections to the inclusive branching ratio BR(6 — > 57), in the case of the MSSM 
with universal boundary conditions, arise mainly from chargino loops and have 
the same sign with the parameter /i. Consequently, these corrections interfere 
constructively with the contribution from the standard model (SM) including 
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4.5 GeV Prediction \i<D 



3.41 GeV tree ievel vaiue ^<0 

3.17 GeV exp. bound 
3.13 tree level value |i>0 



2.3 GeV Prediction n>0 
Exp. bound 2.17 GeV 



Table I. Values for the quark bottom mass. 



an extra electroweak Higgs doublet. However, this contribution is already bigger 
than the experimental upper bound on BR(6 57) for not too large values of the 
CP-odd Higgs boson mass niA- As a result, in the present context with Yukawa 
unification and hence large tan /3, a lower bound on uia is obtained for ^ > 0. 
On the contrary, for < 0, the SUSY corrections to BR(6 57) interfere 
destructively with the SM plus extra Higgs doublet contribution yielding, in 
most cases, no restrictions on the parameters. The results corresponding to the 
parameters given in fig.l are shown in fig. 2. In the case of ;U > 0, BR(6 — > ,37) 
decreases as the splitting between irif^ and increases. 



4. LSP relic abundance and bino stau coannihilations 

The cosmological constraint on the parameter space results from the requirement 
that the relic abundance 0,lsp of the lightest super symmetric particle (LSP) 
in the universe does not exceed the upper limit on the cold dark matter (CDM) 
abundance implied by cosmological considerations {^^lsp is the present energy 
density of the LSPs over the critical energy density of the universe and h is the 
present value of the Hubble constant in units of 100 km sec~^ Mpc"^). Taking 
both the currently available cosmological models with zero/nonzero cosmological 
constant, which provide the best fits to all the data, as equally plausible alter- 
natives for the composition of the energy density of the universe and accounting 
for the observational uncertainties, we obtain the restriction ^lsp 0.22 
(see Refs.||3[)]. Assuming that all the CDM in the universe is composed of LSPs, 
we further get Qlsp h'^ ~ 0.09. 

The cosmological relic density of the lightest neutralino x (almost pure B) 
in MSSM with Yukawa unification increases to unacceptably high values as 
becomes larger. Low values of rriy- are obtained when the NLSP (f2) is almost 
degenerate with x. Under these circumstances, coannihilation of x with f2 and 
f| is of crucial importance reducing further the x relic density by a significant 
amount. The important role of coannihilation of the LSP with sparticles carrying 
masses close to its mass in the calculation of the LSP relic density has been 
pointed out by many authors (see e.g., Refsjl2|, 13, 11)]. Here, we will use 
the method described by Griest and Seckel |11|. Note that our analysis can be 
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Figure 2. The central value of the SUSY inclusive BR(6 sj) as 
function of niA for both signs of fi, Aq = and nif.^ = m^. The contri- 
butions from the SM and the SM plus charged Higgs boson (SM+Higgs) 
as well as the experimental bounds on BR(6 57), 2 x 10^^ and 
4.5 X 10~^, are also indicated. 



readily applied to any MSSM scheme where the LSP and NLSP are the bine 
and stau respectively. 

The relic abundance of the LSP at the present cosmic time can be calculated 
from the equation: 

, 2 1.07 X 10^ GeV"^ 



h 



with 



1/2,^ -1 . 

roo 

aeff = XF {aeffv)x~^dx . 
Jxp 



(7) 

(8) 



Here Mp = 1.22 x 10^^ GeV is the Planck scale, g^, Ri 81 is the effective number 
of massless degrees of freedom at freeze-out and xp = m^/Tp, with Tp being 
the freeze-out photon temperature. 
In our case, a^ff takes the form 



For Tj, we use the nonrelativistic approximation 



ri{x) = 



9eff 



(9) 

(10) 
(11) 



Here gi = 2, 1, 1 (i = X) '^2^ ^2 ) t^*^ number of degrees of freedom of the particle 
species i with mass and x = m^/T with T being the photon temperature. 
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Figure 3. The LSP relic abundance 0,lsp as function of tua in 
the hmiting case m^^ = niy. and for fj, > 0, Aq = 0. The sohd hue 
includes coannihilation of f2 and x, while the dashed line is obtained 
by only considering the LSP annihilation processes. These results are 
affected very little by changing the sign of //. The limiting lines at 
^LSP = 0.09 and 0.22 are also included. 



The freeze-out temperatures which we obtain here are of the order of niy-/ 25 
and, thus, our nonrelativistic approximation (see Eq.(p^) is justified. Under 
these circumstances, the quantities aijV are well approximated by their Taylor 
expansion up to second order in the 'relative velocity', 

aijV = aij + bijv^ . (12) 

The thermally averaged cross sections are then easily calculated 



^3/2 

2^ Jo 



{aijv){x) = / dvv'^{aijv)e 1'^ = + Qhij/x . (13) 



The contribution of the a^'s to the cross section is more important than 
the than the ftj^'s since its contribution is suppressed by a factor Q/xp ~ -2. The 
annihilation cross section is suppressed respect the coannihilation channels due 



to the fact that a^^ is suppressed due to Fermi statistics [10|. This suppression 
is not present in the coannihilation channels, however its contribution to CTeff 
is attenuated by the exponential in rf^ in eq. (10). Therefore coannihilations 
have an important effect in decreasing the ^lsp when the LSP and the NLSP 
are very close in mass. In our case we find this effect negligible when the mass 
splitting x~f is greater than approximately 20%. The complete list of Feynman 
diagrams and the expressions for the a^j's appropriate for large tan/3 are given 
in Ref. §. 

The result of including coannihilation of f2 and x ™- the computation of the 
LSP relic abundance is clearly shown in fig. 3. The effect of increase the mass 
splitting between f2 and x will result in larger values for ^lsp- 
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5. Conclusions 

We have shown that the condition of "asymptotic" Yukawa unification on the 
MSSM results in a significant constraint on the free parameter space of the 
model. 

Constraints from and 6 — > 57 can be simultaneously satisfied for /i > 
and relatively large values of the SUSY parameters ~ 385GeV, ~ 
695GeV, mo ~ 780GeV and M1/2 ~ 1.5TeV. The constraint derived from 
cosmological limits on LSP relic abundance, is satisfied only when bino-stau 
coannihilations are relevant. If we superpose figs. 1, 2, 3 we can observe that 
the previous conditions are satisfied in a narrow band of parameter space. 

The hight values of in the parameter space allowed by the constraints 
described above will make difficult its direct detection. However the large tan /? 
scheme can provide interesting predictions ||T^ if one relaxes the strict unification 
condition imposed in the present study. 
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